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The thermal properties of hydrodynamical simulations of galaxy clusters are usually com- 
pared to observations by relying on the emission- weighted temperature Tow, instead of on 
the spectroscopic X-ray temperature Tspec, which is obtained by actual observational data. 
Here we show that, if the intra-cluster medium is thermally complex, Tow fails at reproducing 
Tspec- We propose a new formula, the spectroscopic- like temperature, T^i, which approximates 
Tspec better than a few per cent. By analyzing a set of hydrodynamical simulations of galaxy 
clusters, we also find that Tsi is lower than Tow by 20-30 per cent. As a consequence, the nor- 
malization of the M-Tsi relation from the simulations is larger than the observed one by about 
50 per cent. If masses in simulated clusters are estimated by following the same assumptions 
of hydrostatic equilibrium and /3-model gas density profile, as often done for observed clusters, 
then the M-T relation decreases by about 40 per cent, and significantly reduces its scatter. 
Based on this result, we conclude that using the observed M-T relation to infer the amplitude 
of the power spectrum from the X-ray temperature function could bias low erg by 10-20 per 
cent. This may alleviate the tension between the value of as inferred from the cluster number 
density and those from cosmic microwave background and large scale structure. 



1 Introduction 

Clusters of galaxies represent a powerful tool for investigating cosmological parameters like the 
normalization of the primordial power spectrum, usually expressed in terms of dg (e.g. Rosati 
et al. 2002). The fundamental quantity entering theoretical models is the cluster mass, M, but 
its high-precision measurement is still precluded by the presence of different systematic effects. 

One way to derive M comes indirectly from the X-ray temperature T, which is related to the 
cluster mass. In fact simple theoretical arguments suggest for virialized systems the existence 



Figure 1: Temperature profiles of tlie simulated galaxy cluster. Filled circles indicate the spectroscopic tem- 
perature profile from Gardini et al. (2003); error bars are at 68 per cent confidence level for one interesting 
parameter. Solid and dotted histogram refer to the mean emission-weighted and spectroscopic-like temperature 
profiles directly derived from the simulation for the same annulus regions. 



a M-T scaling law. The M-T relation is one of the key ingredients in the recipes to extract 
cosmological parameters from the cluster temperature function (XTF) which also requires a 
proper comparison between simulated and observed data. Although simple in principle, this 
comparison is not always trivial. Projection effects and/or instrumental factors introduce sys- 
tematic biases that complicate the correct interpretation of the data. To overcome this difficulty 
we built the X-ray MAp Simulator (X-MAS, Gardini et al. 2004), a software package devoted to 
simulate X-ray observations of galaxy clusters obtained from hydro-N-body simulations. Using 
X-MAS we showed that, if the cluster is highly thermally inhomogeneous, then the projected 
spectroscopic temperature Tspec, obtained from the fit of the observed spectrum, is significantly 
lower than the emission- weighted temperature Te„, which may be directly inferred from the 
simulations by weighting each particle temperature by its emission measure. Depending on the 
thermal complexity of the clusters, this discrepancy may affect the normalization of the M-T 
relation and, consequently, the amplitude of the power spectrum cjg derived from the XTF. 

In Mazzotta et al. (2004) we studied in detail the problem of performing a proper compar- 
ison between temperatures obtained from X-ray observations and simulations. We found that 
although in principle this comparison requires the actual simulation of the spectral properties of 
the clusters, when the cluster temperature is higher than T > 2 — 3 keV, it is possible to define 
a new "spectroscopic-like temperature" function Tgi, that approximates Tgpec to better that few 
percent. In this paper we highlight the main aspects connected with the derivation of Tgi and, 
thanks to a set of hydrodynamical simulations of galaxy clusters, we discuss the implications for 
the M-T relation and for the determination of ag . 

2 Spectroscopic-like temperature 

From a purely analytic point of view the spectroscopic temperature Tgpec is not a well defined 
quantity. In fact, it results from the fit of a single-temperature thermal model to a multi- 
temperature source spectrum; however, since the former cannot accurately reproduce the spec- 
tral properties of the latter, it follows that Tgpec cannot be unequivocally identified. Despite of 
this, as observed spectra are affected by a number of factors that distort and confuse some of 
their properties, in some circumstances multi-temperature thermal source spectra may appear 
statistically indistinguishable from single-temperature ones. In Mazzotta et al. (2004) we dis- 
cuss this aspect by focusing our attention on observations of clusters of galaxies made using the 
CCD detectors of Chandra and XMM- iVewton. From this study we find two very important 
results that we summarize as: 



i) given a multi-temperature source spectrum, if the lowest dominant temperature component 
has Ti > 2 — 3 keV, then a fit made with a single-temperature thermal model is statistically 
acceptable regardless of the actual spread of the temperature distribution of the source; 

ii) the emission- weighted temperature Tgw does not properly estimate the spectroscopic temper- 
ature Tspec as it tends to overestimate it. The larger the spread in temperature of the thermal 
components of the observed spectra, the larger the discrepancy. 

To partially solve this problem we propose to use the spectroscopic-like temperature function 
Tgi instead of Tew. The idea behind the analytic derivation of T^i is quite simple and is itemized 
below. 



Given a multi-temperature source spectrum we are interested in identifying the temper- 
ature Tgi whose spectrum is closest to the observed one. If for simplicity we consider a 
source spectrum composed by two thermal models with densities ni, n2 and temperatures 
Ti, T2 respectively, requiring matching spectra means that: 



nlaZ, Ti)^ exp(-^) + njCiZ, T2, )^ exp(-^; 
«^C(^,T..)^exp(-^), 



(1) 



where A is an arbitrary normalization constant and C{Z, T) is a parameterization function 
accounting for the total Gaunt factor and partly for the line emission. We approximate 
C{Z,T,E)cx{T/keYr. 

Supported by the fact that both Chandra and X.MM-Newton are more sensitive to the 
soft region of the X-ray spectrum, we can expand the previous equation in Taylor series 
to the first order of E/kT and obtain: 



where we generalized to an arbitrary number of thermal components for the source spec 



Tsi = 7^^— — 372 ' (2) 



trum. 



• The power index a needs to be calibrated to the specific detector. This can be done by 
simulating a set of two-temperature source spectra and minimizing the value of the mean 
percentile variation of Ts\ with respect to Tspec- In Mazzotta et al. (2004) we show that 
a = 3/4 represent a good choice for each CCD detector of Chandra and XMM.-Newton. 

It is important to say that we tested Tgi against Tgpec, and Tew by using the hydro-N- 
body simulation of a thermally complex (major merger) cluster of galaxies and by applying 
X-MAS. We compared the projected temperature maps and profiles derived by using different 
temperature definitions. We find that Tgi approximate Tspec to better than few percent, while 
Tew may overestimate Tgpec by factors as large as 2 (see Mazzotta et al. 2004 for details). Just 
as an example in Fig. ^ we show the azimuthally averaged temperature profiles of this cluster 
as obtained by adopting the three different estimators. It is worth noticing that while Tgy is 
consistent within the errors with Tgpec, Tew is always higher. 



3 Cosmological implications 

In this section we discuss the consequences of the discrepancy between Tgpec and Tew for the 
M-T relation and for the estimate of erg from the XTF. To do that we present the analysis of 
a sample of hydrodynamical simulations of galaxy clusters for which we use the temperature 



function Tsi defined in the previous paragraph. The results presented here are extracted from 
Rasia et al. (2005) to which we refer for more details. 

We created the sample of simulated galaxy clusters by merging two different sets: 



• Set#l: 95 temperature-selected clusters (with Tew > 2 keV), extracted from the large- 
scale cosmological hydro-N-body simulation of a "concordance" ACDM model (llom = 0.3, 
OoA = 0.7, Qob = 0.019 /^-^ h = 0.7, ag = 0.8), presented in Borgani et al. (2004, B04) 

• to account for the high temperature (T > 7 keV) clusters not present in Set#l because of 
the limited box size of the previous simulation, we added a Set#2 composed by clusters 
having Mvir > 10^^h~^ Mq belonging to a different set of high-resolution hydro-N-body 
simulations (Dolag et al. 2005 in preparation). 

Besides gravity and hydrodynamics, both simulations includes the treatment of radiative cooling, 
the effect of a uniform time-dependent UV background, a sub-resolution model for star forma- 
tion from a multiphase interstellar medium, as well as galactic winds powered by SN explosions 
(Springel & Hernquist 2003). 

For each cluster in our sample we calculate both Tew and Tgi. In Fig. [3 we plot the values 
of Tsi versus Tew. The distribution of the points in Fig. [2 suggests that the two estimators can 
be roughly related by a linear relation. By applying a fit; we find T^i = (0.70 ± 0.01)Tew + 
(0.29 ± 0.05), shown as dotted line in Fig. HI 

The existence of a systematic trend between the two different estimates of the cluster tem- 
perature has interesting consequences on the M-T scaling relation. Using Set^T^l, Borgani et al. 
(2004) showed that the normalization (Mq) of the simulated M-Tew relation is 20% higher than 
the observed one. Since Tew tends to overestimate Tgi by a fair amount, we expect the actual 
M-Ts\ relation from simulations to be even more discrepant with respect to observations. In 
fact, we find that Mq is higher than the observed one by about 50 per cent. However, in order 
to compare the M-T obtained in observations and simulations, it is also necessary to derive 
the simulated clusters mass under the same assumption of the X-ray analysis, i.e. by applying 
the condition of hydrostatic equilibrium to a spherical gas distribution described by a /3-model 
(Cavaliere & Fusco-Femiano 1976), with the equation of state having a polytropic form. In this 
way, the total self-gravitating mass within the radius r is given by (Finoguenov et al. 2001; 
Ettori et al. 2002), 

M« r) = Lll X 10"ft.7^p4^^ /.-'Me. (3) 

Here T(r) is the temperature at r, /?fit is the fitted slope of the gas density profile, x = r/rc is the 
radial coordinate in units of the core radius Tc, and 7 is the effective polytropic index obtained 
from temperature and gas profiles. Fitting the simulation results gives log(Mo//i~^M0) = 
13.41 lb 0.13 when fixing a = 1.5. 

As already mentioned, the M~T relation is one of the key ingredients in the recipes to 
extract cosmological parameters from the cluster XTF: under the reasonable assumption that 
the scaling relations do not significantly change with (Tg, a larger Mq implies a larger mass for a 
fixed temperature and, therefore, a higher normalization of the power spectrum (e.g., Borgani et 
al. 2001; Seljak 2002; Pierpaoh et al. 2003; Henry 2004). Huterer & White (2002) suggested an 
approximation for the scaling of (Tg with Mq, involving the matter density parameter: r2o;^(Tg oc 
Mg '^^. Based on this relation and on the fact that a number of analyses of hydrodynamical 
simulations (e.g., Muanwong et al. 2002; B04; Rasia et al. 2004; Kay et al. 2004) have shown 
that EqEl underestimates the actual cluster mass by about 20%, we expect that the value of dg 
is underestimated by a similar amount, i.e. about 20%. 
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Figure 2: The comparison between TLw and Tsi. Circles 
and triangles refer to objects from Set^^l and Set#2, 
respectively. The dotted line is the best fit, while the 
dashed line corresponds to Tgi — Tew 
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Figure 3: The cumulative temperature function N{> 
T), from Set#l. Filled squares and dashed line refer 
to the results for Tew, while filled triangles and solid 
line represent T^i- Error bars correspond to Poissonian 
uncertainties. For reference, open circles show the local 
temperature function measured by Ikebe et al. (2002). 



In FigEl we show the effect on the XTF of using Tgi, instead of Te„. We compare here 
the cumulative temperature function A^(> T) from the simulation (Set#l only) to the data by 
Ikebe et al. (2002). The squared points represent the emission- weighted temperature; being 
Tew in good agreement with the data, one could argue that the as used in the simulation well 
represents the observed universe. This would imply as = 0.8. However, as already shown, 
the emission-weighted temperature is typically an underestimate of the observed temperature, 
while the spectroscopic-like temperature better takes into account the observational procedure 
of temperature measurement. The comparison between Tgi and data shows that as = 0.8 is 
excluded at 3a level, while favoring a larger power spectrum normalization, at least for the ICM 
physics included in our simulation. Computing the mass function by Sheth & Tormen (1999) 
for our cosmology at the mass corresponding to a 2 keV cluster, we find that as needs to be 
increased to about 0.9 to increase the number density of clusters above this mass limit by the 
required 50-60 per cent. 



4 Conclusions 

In this section we simply summarize the main results presented in this paper. 

1. The projected spectroscopic temperature Tgpec of thermally complex clusters obtained 
from X-ray observations is always lower than the emission- weighed temperature Tew, which 
is instead widely used in the analysis of numerical simulations. We show that this temper- 
ature bias is mainly related to the fact that the emission-weighted temperature does not 
reflect the actual spectral properties of the observed source. 

2. A proper comparison between simulations and observations needs the actual simulations 
of spectral properties of the simulated clusters. Nevertheless, if the cluster temperatures is 
> 2 — 3 keV it is possible to define a temperature function, that we call spectroscopic-like 
temperature Tgi, which approximate Tgpec to better than few per cent. 

3. By analyzing a sample of hydrodynamical simulations of galaxy clusters, we find that Tgi 
(and, thus Tgpec) is lower than Tew by 20-30 per cent. We obtain a linear fit approximat- 
ing the relation between the two different temperature estimators. As previous studies 
made using Tew show that the discrepancy in the M~T relation between simulations and 



observations is about 20 per cent, it is clear that the use of Ts\ increases this discrepancy 
to ~ 50 per cent. Nevertheless, if we assume hydrostatic equilibrium for the gas density 
distribution described by a /3-model with a polytropic equation of state, we know that 
masses are underestimated on average by ^ 40 per cent. Although this goes in the direc- 
tion of substantially reducing the discrepancy with observational data, this is not sufficient 
to cancel it. 

4. The bias in the M-T relation propagates into a bias in (jg, from the XTF. If such a bias 
is as large as that found in our simulations, the values of cjg obtained by combining the 
local XTF and the observed M-T relation are underestimated by about 15 per cent. 

5. The XTF from the simulation is significantly lower when using Tgi instead of Tew. A 
comparison with the observed XTF indicates that for the "concordance" ACDM model as 
needs to be increased from ~ 0.8 to ~ 0.9. 

To conclude, the results of this study go in the direction of alleviating a possible tension 
between the power-spectrum normalization obtained from the number density of galaxy clusters 
and that arising from the first-year WMAP CMB anisotropies (e.g. Bennett et al. 2003) and 
SDSS galaxy power spectrum (e.g. Tegmark et al. 2004). 
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